multistage pcy algorithm